Integral relations for three-body continuum states with the adiabatic expansion 



OS 

o 
o 

(N 

Or 
<D 



O 
=5 



(N 
> 

in 
O 
(N 

in 
o 
o> 
o 



P. Barletta, 1 C. Romero-Redondo, 2 A. Kievsky, 3 M. Viviani, 3 and E. Garrido 2 
1 Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom 
2 Institute de Estructura de la Materia, CSIC, Serrano 123, E-28006 Madrid, Spain 
3 Istituto Nazionale di Fisica Nucleare, Largo Pontecorvo 3, 56100 Pisa, Italy 

Application of the Hyperspherical Adiabatic expansion to describe three-body scattering states 
suffers the problem of a very slow convergence. Contrary to what happens for bound states, a huge 
number of hyperradial equations has to be solved, and even if done, the extraction of the scattering 
amplitude is problematic. In this paper we show how to obtain accurate scattering phase shifts 
using the Hyperspherical Adiabatic expansion. To this aim two integral relations, derived from the 
Kohn Variational Principle, are used. The convergence of this procedure is as fast as for bound 
states. 
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Introduction. — Few-body collisions involving either 
nuclei, atoms, or molecules are frequently investigated. 
To this aim different methods are at present available 
depending on the interaction under study. In nuclear 
physics, collisions involving three or four nucleons have 
been extensively studied within the Faddeev method or 
the Hyperspherical Harmonic (HH) method [l], 0, 
These two methods show sufficient flexibility to treat 
the complexities of the nucleon-nucleon interaction. A 
different problem arises when the interaction presents a 
hard core, as in the case of the atom-atom interaction, 
or in systems with A > 4. In the first case, the Faddeev 
equations have been extended to deal with a hard core 
repulsion 0] whereas the Hyperspherical Adiabatic (HA) 
expansion method proved to be a very efficient tool 
In nuclear systems with A > 4 tentatives to describe 
scattering states have recently appeared 0, 0] ■ 

Here we are interested in describing a 1+2 collision 
using the HA expansion method. For bound states the 
convergence of the HA expansion has proved to be very 
fast. However the convergence of the expansion slows 
down significantly in the case of low energy scattering 
states Q. On the other hand this method is extensively 
used to describe few-atom systems in the ultracold regime 
(see Refs. d, [l(| and references therein) and, in particu- 
lar, atom-dimer collisions. Therefore a detailed study of 
its convergence properties is timely. 

In this letter it is shown for the first time how the HA 
expansion method can be used to describe elastic scat- 
tering with a pattern of convergence similar to a bound 
state calculation. This is achieved in a simple but very 
general way in which a second order estimate of the phase 
shift is extracted from the wave function using two inte- 
gral relations derived from the Kohn Variational Prin- 
ciple (KVP) [ll|. The number of HA terms needed to 
obtain completely stable results depends very little on 
the structure of the potential, exactly as for bound state 
calculations. The integral relations are governed by the 
wave function in the interaction region. Therefore the 
stability of the results with a low number of HA basis el- 
ements is a clear indication that inclusion of more terms 
in the expansion only modifies the wave function outside 



the interaction range. 

As derived from the KVP, the integral relations are 
general and their application is not limited to three par- 
ticles. They can be applied to an >l-body system in which 
the scattering wave function is known in the interaction 
region. Example of applications are given below. 

Continuum states in the Hyperspherical Adiabatic Ex- 
pansion method. — The details of the HA method can 
be found in 0, Q . For simplicity, here we restrict our- 
selves to three equal mass particles with total angular 
momentum L = and with only s-waves involved. 

From the Jacobi coordinates Xi = (rj — rfc)/v2 and 
Y i = (rj + r& — 2r i )/v / 6, one defines the hyperspher- 
ical variables, [p,fl*] = [p, <f>i,Pi], with //, = X, ■ Y i} 
Xi = p cos <f>i, and Yi = psin^, where {i, j, k} is a cyclic 
permutation of {1,2,3}, and {r^} are the coordinates 
of the three particles. In hyperspherical coordinates the 
Hamiltonian operator Ti. takes the form: 
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.G 2 + V( P ,n)=-—T p +Hn, (1) 
2m 



where T p is the hyperradial operator, G 2 is the grand- 
angular operator, V(p,fl) = ^iV(Xi) is the potential 
energy, and m is set equal to the mass of the particles. 
The wave function '5 for a specific bound or continuum 
state is expanded as: 



STTl 



(2) 



where S, T, and n are the total spin, total isospin and 
parity. For simplicity we shall suppress from now on the 
corresponding labels in ^ and The HA basis ele- 

ments are the eigenfunctions of Hq at fixed values 

of p. Their corresponding eigenvalues, U v (p), are the 
adiabatic potentials, which enter in the coupled set of 
differential equations (see Refs. (H. Hojl 

| -& t p + Mp) - &Q»Ap) - e] u u {p) - 

N A 

[Q^'(p)+Puu'(p)(l + 2^)]u u >(p) = (3) 
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with Na the number of adiabatic channels included in 
the calculation, E the three-body energy, and from which 
the hyperradial functions u v (p) are obtained. At energies 
below the two-body breakup E 2B and p — > oo, the total 
scattering wave function behaves asymptotically as [8j 
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\ST). (4) 



However, as we will show below, even increasing Na as 
much as possible, the computed value of 5 P does not con- 
verge to the expected one. This can be understood from 
the fact that the asymptotic structure of the system can 
be constructed in terms of the functions: 



F ST = £ F ST (i) = £ MXd^P^lST) 



i=l 
3 



i=l 
3 



Gst = £ G ST « = £ MXi)^^-\ST: 



(5) 



where particle i is assumed to hit the bound state made 
by j and k, and where yi = ^Yi is the distance be- 
tween i and the j-k center of mass, and k 2 = lk 2 „. The 
- 1 y o p 

asymptotic configuration in the limit yi — > oo is then: 



Fst + tan S Gst- 



(6) 



The function Gst, representing a regularization of the 
function Gst, introduces a nonlinear parameter 7 to 
eliminate a term proportional to 5(y i ) originated by 
(H — E)Gst, and fy c is the part of the wave function in- 
side the interaction region constructed in terms of some 
parameters (e.g. a linear combination of basis elements). 
It verifies & c — > asymptotically. The other parameter 
in \1/ is tan S . In the present work we have used 

Gst = £ UXi) C °^ Vi \ l ~ e~^)\ST). (9) 
i y k y yi 

The variation of the functional (JSj> with respect to the 
parameters in \& c and with respect to tan S leads to: 

<* c |£|*>=0 ; < G ST \C\^ >= 0. (10) 

These two equations can be interpreted in two different 
ways. In the case in which \& is explicitly separated in 
the three terms Fst, Gst, the above equations are 
used to determine ^ c and the first order estimate of tan 5 
(tan<5 13 ). Accordingly, * is constructed after solving 
these equations. A different case arises when ^ is known 
(for example using the HA expansion) but the separation 
in the three terms is not explicitly known. For this case 
the two equations can be used to define of ^ c and tan S 1 . 

Introducing Eq. (JTDJ) into the functional, the second 
order estimate of tan <5 is obtained: 



When p — * 00, the distance Xi is limited by (\>d and 
the approximate relation k y yi ss k p p holds. However, the 
exact equivalence between k y yi and k p p is not matched 
for any finite value of p and, accordingly, the boundary 
condition of Eq.(d]) is equivalent to the one in © only 
at p m 00 and Na. — > 00. As a consequence S p converges 
extremely slowly to 5 by increasing the number of adia- 
batic states. This situation has reduced the applicability 
of the method. 

Second order integral relations. — From the above 
discussion and observing that in the expansion of the 
functions Fst and Gst in terms of HA basis elements [1] , 
the two terms of Eq. (fj| represent the first term of that 
expansion, respectively, the wave function ^ can be ex- 
pressed asymptotically as: 
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[tan 5] 



(tan<5) ist - < F ST \C\{1/A)^ > (11) 



with A =< ty\£\GsT >■ This is a consequence of the 
general relation A =< ^\C\G S t > - < G S t\C\^ > (ob- 
tained by transforming the Laplacian term in a surface in- 
tegral), the normalization relation < Fst\£\Gst > — < 
Gst I £ I Fst >= 1 and the last equation derived from the 
KVP in Eq. (fT0| . The same relation can be used to ob- 
tain a first order estimate for the coefficient B as 



-I St 

B 1 = 



< F ST |£|* > - < f\C\F S T > 



(12) 



After multiplying Eq. (TlTj) by A one gets that B 2 " = 
B ist - < F S t\£\^ >, which by use of Eq. {T2J) leads to 
a second order integral relation for B and, accordingly, 
a second order estimate for tan S. These results are the 
main conclusions of this paper and can be summarized 
as 



In order to extract the coefficients A and B (tan 8 = 
B/A), we derive from the KVP two integral relations 
accurate up to second order. The KVP states that the 
following functional is stationary: 

[tan<5] 2 " d = tim5- < (1/A)V\C\(1/A)$> > (8) 



with respect to variations of the wave function, where 
C = -^yp(Tl — E). The scattering wave function can be 

schematically written as (1/A)^ = ^c+FsT+tanJ Gst- 



B 2 " d = - 
A 



-<f4^> jtan^=B 2 ">. (13) 
<V\£\Gst> J 1 y 1 



The relations of Eq. (fl3)) are equivalent to the KVP and 
are useful in the cases in which ^ is known but its explicit 
asymptotic form in terms of the functions Fst and Gst 
is not. This is the case, for example, when * is obtained 
from the solution of the HA equations. The integrands 
in the integral relations of Eq. (TT5|) go rapidly to zero as 
p — > 00 since -Fst and Gst are solutions of C in that 



3 



limit. Therefore an accurate knowledge of ^ outside the 
range of interaction is not needed. In the present case, 
the explicit form of the integrals in Eq. (fT3|) are: 

B 2 " d =-C J dpfPdnVfatoMXiXFsrW + FsTik)] 

A = C J dpp 5 dn*( P ,Sl)V{Xi)[GsTti) +G S T(k)} +7 7 

where C = 2\/3m/h 2 and J 7 is a (short-range) integral 
including all terms depending on 7. Let us note that the 
last integral is largely independent of 7 provided that the 
regularization is performed inside the interaction region 
and $ tends to the exact wave function. The dependence 
of tan 5 o n 7 is studie d below. We have found that values 
of 7 w y / rn\E^\/W are a convenient choice. 

The integral relations, as given above, are a gen- 
eralization of the relation used in the two-body case 
[l2l ]. However, an attempt to identify the single term 

B 2 — tan as a corrected phase-shift fails as we show 
below. The validity of Eq. (TlB"]) is not limited to the three- 
body case and to wave functions obtained using the HA 
method. It can be applied to any wave function '5 which 
verifies {TL — E)^ — in the interaction region without 
any explicit indication of its asymptotic behaviour. An 
example is represented by the solution of (H — E)^f = 
in a box in which ^ is set to zero at some distance. Using 
Eq. (fTB")) a second order estimate of a phase shift can be 
obtained studying its convergence in terms of the dimen- 
sion of the box. 

Results. — As a first application we consider a three- 
body system of identical spinless bosons interacting 
through a central, s-wave, gaussian potential V(r) = 
Vq exp — (r/ro) 2 , with Vq — —51.5 MeV and r = 1.6 
fm. Though this potential is unrealistic, it will serve 
to the purpose of testing the method due to the very 
extended dimer wave function (E 2B = 0.397743 MeV 
and h 2 /m = 41.4696 MeV fm 2 ). Such three-body sys- 
tem has two L=0 bound states with separation energies 
E$ = -9.7574 MeV and E^, = -0.4816 MeV, respec- 
tively. In Table [T] we show the convergence of these two 
states in terms of the HH and HA expansions. From the 
table we observe the fast convergence of the HA even in 
the case of the shallow state E^, . We can conclude that 
10 HA basis states are sufficient to describe simultane- 
ously both bound states. 

We now show results for the L=Q phase shift at E = 
-0.1 MeV. Eq.Q has been solved up to p=500 fm with 
the boundary condition of Eq. ([J} for increasing values of 
Na- Due to the large extension of the P„„/ and Q vv > cou- 
pling terms, the higher N a the larger the value of p at 
which the asymptotic form of u\ is verified. For example, 
using only one HA basis state (only one hyperradial equa- 
tion has to be solved), U\ reaches its asymptotic form at 
p sa 100 fm. When 40 HA terms are used, this happens 
beyond 500 fm. The results for S p , <5e=arctan(i? 2 " ), 

S 2 = arctan(£> 2 /A) (for different values of 7), and A, 
are given in Table HT1 up to 40 HA basis functions. 
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HH 


HA 


N 


HH 


HA 


1 


-9.2062 


-9.7347 


1 




-0.4781 


2 


-9.5810 


-9.7552 


4 




-0.4815 


3 


-9.7247 


-9.7573 


10 


-0.2323 


-0.4816 


4 


-9.7424 


-9.7574 


30 


-0.4635 


-0.4816 


6 


-9.7558 


-9.7574 


50 


-0.4790 


-0.4816 


8 


-9.7571 


-9.7574 


70 


-0.4811 


-0.4816 


10 


-9.7574 


-9.7574 


100 


-0.4815 


-0.4816 



TABLE I: Convergence of the bound state energies (in MeV) 
as a function of the number TV of HH and HA basis functions. 
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B 




A 




5 2 ' 
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arctan(73 2 " d /A) 










7 = 


-- 0.1 


7 = 


0.25 


7 = 0.5 


7 = 


= 1.0 


1 


65 


23 


67 


85 





607 


75 


895 


76 


212 


75.479 


74 


727 


4 


71 


65 


71 


89 





909 


73 


385 


73 


446 


73.450 


73 


429 


8 


72 


32 


72 


42 





951 


73 


222 


73 


237 


73.245 


73 


259 


12 


72 


56 


72 


62 





965 


73 


190 


73 


194 


73.196 


73 


203 


16 


72 


67 


72 


71 





971 


73 


182 


73 


183 


73.185 


73 


189 


20 


72 


72 


72 


76 





974 


73 


180 


73 


180 


73.182 


73 


186 


24 


72 


75 


72 


79 





976 


73 


179 


73 


179 


73.180 


73 


184 


28 


72 


77 


72 


81 





977 


73 


179 


73 


179 


73.180 


73 


184 


32 


72 


78 


72 


82 





978 


73 


179 


73 


179 


73.180 


73 


184 


36 


72 


79 


72 


82 





978 


73 


179 


73 


179 


73.180 


73 


183 


40 


72 


79 


72 


83 





978 


73 


179 


73 


179 


73.180 


73 


183 



TABLE II: Patterns of convergence for 5 P , 5b = arctanB, 

S 2 (in degrees) and A, in terms of Na for E=— 0.1 MeV. 
The values of A have been calculated using 7 = 0.25. 



We observe that S p and 5b converge very slowly to a 
value that, by extrapolation, can be estimated in the in- 
terval 72.8° — 72.9° as Na — > 00. This is at variance with 
the value S 2 = 73.18°, which shows a rate of conver- 
gence extremely fast and a large stability with 7, as Na 
increases. The calculation of S 2 requires the knowledge 
of the radial functions up to values of p not larger than 
70 — 80 fm, for which a relatively small number of HA 
terms is enough. Conversely, 5 p and 5b would converge 
to the correct phase shift only after imposing the bound- 
ary condition to the wave function at p = 00, for which in 
principle infinitely many HA basis terms are needed. For 
comparison, a converged value of 73.180° is obtained for 
the phase-shift using the HH expansion method with 120 
basis elements. This result is in complete agreement with 
the one obtained with Eq. (fTB)) . The different patterns of 
convergence can be clearly seen on Fig[T] 

Ref. Q reports calculations at three different energies 
using the MT-III potential [13j, which has a yukawian 
repulsion at short distances, in the S = 3/2, T = state. 
When the HH expansion is used, more than 120 basis 
states have to be included to reach convergence in the 
phase shifts. In Tabic UTTI we show the corresponding re- 
sults when using Eq. (| 1 3|) . In the last row the results from 
Ref. [8j using the HH expansion are given for compari- 
son. The results obtained with the integral relations are 
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FIG. 1: Convergence of S p , 5b, and S 2 as a function of 
the number Na of HA basis functions. The converged value 
obtained from the HH expansion is shown for comparison. 



order. Their ratio, the phase shift, converges in terms 
of the HA basis elements as fast as the binding energy 
in a bound state calculation. The fast convergence has 
been shown for different types of interactions. We would 
like to stress the general validity of Eq. (p~3|) . Its applica- 
tion will be very useful in the case in which 'J is known 
in the interaction region but the exact construction of 
its asymptotic form is difficult. The HA expansion has 
been applied in Refs. [l(| EH to compute phase-shifts in 
a 1+2 and a 2+2 helium atom collisions. Accordingly. 
Eq. (fT3| can be used directly to obtain a second order 
estimate of the phase-shifts. In Ref. Q, n — a scatter- 
ing has been studied using Quantum Monte Carlo tech- 
niques. The wave function of the system was obtained 
solving (Ti — E)*i> = in a box. The knowledge of W 
in the interaction region allows for a direct application 
of Eq. (|13[) also in this case. To be noticed that in the 
case in which more than one elastic channel is open, the 
coefficients A and B of Eq. (fl~3|) correspond to matrices 





0.2 MeV 


1.0 MeV 2.0 MeV 


4 


-28.277 


-55.875 


-71.507 


8 


-28.290 


-55.865 


-71.475 


12 


-28.293 


-55.864 


-71.473 


16 


-28.294 


-55.863 


-71.473 


20 


-28.294 


-55.863 


-71.473 


HH 


-28.294 


-55.863 


-71.474 



TABLE III: Convergence of S 2 " = arctan(£ 2 " /A) (in de- 
grees) at three incident energies with the MT-III potential. 



in complete agreement with those obtained using the HH 
expansion and show a very fast pattern of convergence. 

Conclusions. — We have derived two integral rela- 
tions from the KVP which are accurate up to second 



Bf = - < 9^ > 
™ = <%\C\G j > 



(14) 



R 2 = A~ 1 B 2 " 

A_ 

with R 2 " the second order estimate of the scattering ma- 
trix whose eigenvalues are the phase shifts and the inde- 
ces indicate the different asymptotic configurations 
accessible at the specific energy under consideration. Fi- 
nally we would like to mention the possibility of using 
Eq.(13) to describe a 1+2 elastic collision with charged 
particles using a screened Coulomb potential and free 
asymptotic conditions. 
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